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Recently, Kozik, Ferrero and Georges have discovered numerically that for a family of fundamental 
models of interacting fermions, the self-energy 5][G] is a multi-valued functional of the fully dressed 
single-particle propagator G, and that the skeleton diagrammatic series Sboid[G] converges to the 
wrong branch above a critical interaction strength. We consider the zero space-time dimensional 
case, where the same mathematical phenomena appear from elementary algebra. We also find a 
similar phenomenology for the fully bold formalism built on fully dressed single-particle propagator 
and pair propagator. 


In quantum many-body physics, an important role is played by the self-consistent field-theoretical formal¬ 
ism, where the self-energy S is expressed in terms of the exact propagator G (see, e.g., [1] and Refs, therein). 
In a recent article, Kozik, Ferrero and Georges numerically discovered mathematical difficulties with this 
formalism [2]. They studied not only the Hubbard model in two space dimensions, but also simpler models 
—the Hubbard atom and the Anderson impurity model— for which there is no spatial coordinate. Here we 
consider an even simpler toy-model for which there is no imaginary-time coordinate either. This idea was 
also followed in the very recent article [3]. For fermionic many-body problems, the partition function can be 
written as a functional integral over Grassmann fields in {dpi) space-time dimensions {d spatial coordinates 
and one imaginary-time coordinate) [4]. Accordingly, we consider the zero space-time dimensional model 
defined by a “partition function” given by a simple Grassmann integral. 
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with the action 


<J 


( 2 ) 


and a corresponding propagator 


G — {(pfjcpfj) 
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Here a G {t, 1} is the spin index, while /i and U are dimensionless parameters that play the roles of chemical 
potential and interaction strength. Diagrammatically, the Feynman rules for the present toy-model are 
analogous to the ones of the physical (d + 1) dimensional models, with the simplification that space-time 
variables are absent and that the propagators are constants. 

In this exactly solvable toy-model, we observe a similar phenomenology than the one found by Kozik et 
al. in non-zero space-time dimensions. More precisely, restricting to U < 0, we find that: 

• The mapping Go G(Go, U) is two-to-one and hence the function G i-G- S(G, U) has two branches. 


• The skeleton series Sboid(G, G), evaluated at the exact G(/i, I/), converges to the correct branch for 
\U\ < /i^, and to the wrong branch for \U\ > jj?. 


This can be derived very directly from the above definitions, 
yields 


Expanding the exponentials in 


Eqs. ^ 


Z{ii, U) = II^-U 


( 4 ) 



G{fi,U) = 




n'^-u 


The propagator for = 0 is 


Go(m) = -• 


(5) 
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The self-energy S, defined as usual by the Dyson equation G ^ = Gq — S, reads 

li 

We note that for [/ > 0, an obvious pathology appears in this model around U = namely, Z changes 
sign, and G diverges. Therefor e we restrict to U <0. 

Eliminating fi between Eqs. (5|6) gives 


G(Go,t/) = ^ 


Go 


UGl 


(8) 


The map Go G(Go, U) is two-to-one, because the Go’s that correspond to a given G are the solutions of 
the second order equation 


GGGg + Go-G = 0, 


which has the two solutions 


Gf\G,U) = 


-l±\/l+4f/G2 
2UG ■ 


These solutions are real provided (G, U) belongs to the physical manifold {{G{ij,, U), U)}; indeed, 

4\U\G{p,Uf <1. 

The corresponding self-energies (given by the Dyson equation) are 


U) = 


-1± V1 + 4C/G2 
2G ■ 


(9) 


( 10 ) 


( 11 ) 


( 12 ) 


The correct self-energy is recovered from T,^‘^\G{ii,U),U) provided one takes the determination 

s = sign(/i2 - \U\). (13) 

We turn to a discussion of the skeleton diagrammatic series Ilboid(G, G) for the self-energy E in terms 
of fully dressed propagator G and bare vertex U. We find that Eboid(G, G) is the U^O Taylor series of 
E^+^(G, G). Before deriving this, we note that obviously, Eboid(G, G) can never be the Taylor series of 
E(-)(G,G), since the former vanishes at G = 0 while the latter does not. Eor the derivation, it is convenient 
to introduce g := y^]U\G and go := so that Eqs.(8 10) simplify to g{go) = go/{l + 


96 


o^\g) = (1 T a/I — 4:g‘^)l{2g). The key point is that go^\g{go)) = 90 ^ where the symbol = means equality 
in the sense of formal power series. This is because the inverse mapping of g{go) is go^\g) for small go and g. 
Let us then denote by Ebare(Go, U) the diagrammatic series for the self-energy in terms of bare propagators 
and vertices. Setting Eboid(G, 17) =: y/\U\ ahold{g) and I]bare(Go, 17) =: vTi'^bare(5o)> a defining property 
of cTboid is that crboid(5(5o)) =crbare(5o)- In the present toy-model, we simply have Sbare(Go, 17) = 17Go, i.e., 

ab.M = - 90 . Hence, (Tboid(<7)= - ie., Ehoid{G, U) = U Gi+\G,U) = S(+)(G,f/). 

Explicitly, expanding the square root in Eq. ([T^ yields 


Sboid(G,f/) = y] 


n=l 


(-l)n-i( 2 n- 2 )! 

T.! (n — 1)! 


Q^n—ljjn 


(14) 
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FIG. 1. The two branches of the self-energy E, multiplied for convenience by G, as a function of the interaction 
strength \U\, for hxed fi. Dotted line: exact self-energy Eexact = dashed line: skeleton series Eboid(G, f/). 

The upper branch corresponds to E^^^(G, f/), the lower branch to U). Here, G stands for the exact G(/i, U). 


It is natural to evaluate the bold series at the exact G{iJi^U). The obtained series always converges, as 
follows from the inequality The convergence is always to E^+^(G(/i, t/), t/), which as we have seen is 

the correct result for \U\ < the wrong one for \U\ > . The convergence speed is slow for \U\ close 

to and gets faster not only in the small \U\ limit, but also in the large \U\ limit. This is qualitatively 
identical to the numerical observations of Kozik et al in non-zero space-time dimensions. We note that the 
series converges even at the critical value \U\ = albeit very slowly (the summand behaving as 1/n^/^ for 
large n); at this point, the boundary of the series’ convergence disc is reached. 

In the Figure we plot the quantity E G, which, for the exact E, is equal to U times the double occupancy 
versus \U\ for fixed p. The picture is qualitatively identical to Fig. 2(a) of Kozik et al. 

Geometrically, the mapping U T^{G,U) can be viewed as single-valued on a two-sheeted Riemann- 
surface with a branch point at —1/(4G^). Let us vary U from 0 to —oc for fixed p. For small |I7|, the point 
U is far away from the branch point and the bold series converges quickly. The result corresponds to the 
correct Riemann-sheet. Upon increasing |[/|, the point U and the branch point — l/(4G^(/i, U)) both move 
leftwards. The point U catches up the branch point when \U\ = p^. For larger |U|, U is again to the right 
of the branch point, and the bold series converges again, but the result corresponds to the wrong sheet. In 
principle, the correct result can be recovered from Eboid(G, U) by analytic continuation along a path where 
U rotates once around the branch point. 

As emphasized by Kozik et a/., since the self-energy is the functional derivative of the Luttinger-Ward 
functional 4>[G,U] with respect to G, and since E[G, U] is multivalued, 4> must also be multivalued. This 
can also be seen explicitly in the present model. The Luttiger-Ward functional (which is actually a function 
in the present model) can be constructed following the usual procedure (see, e.g.^ [T]). Starting from the free 
energy F(/i, U) = — lnZ(/i, U), and noting that 


dF 

dp 


(/r,U) = -2G, 


the Baym-Kadanoff functional is defined by Legendre transformation: 


n{G,U) = F{p,U)F 2 pG 


with /i(G, U) such that Eq. (15) holds. The Luttinger-Ward functional is then defined by 


4>(G,U) = O(G,U)-U(G,0). 


(15) 


(16) 


(17) 


This leads to the expression 

^±{G,U) = -ln|S(=F)(G,[/)| -ln|G| -2GI](=F)(G,f/) -2. (18) 
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There are two branches because Eq. (15) has two solutions. Accordingly, the mapping /i i-^ E(/i, U) is neither 
convex nor concave. Finally one can check thal[^ 




{G,U) = 2^^^\G,U). 


(19) 


We point out that in the present zero space-time dimensional model, both branches (G, U) are physical 
in the sense that they are the non-interacting propagator for certain parameters of the model. This is not 
the case for the Hubbard atom and the Hubbard model [2] . More generally, the absence of imaginary-time 
coordinate consitutes a drastic simplification, and while we have observed similar phenomena than in I2],it 
remains an open question to which extent the underlying mechanisms are similar. 

Finally we briefly treat the fully bold formalism built not only on the fully dressed G but also on the fully 
dressed pair propagator T. This formalism was used, e.g., in Refs. One defines 

T = U -U^ (20) 

or diagrammatically 

( 21 ) 



The pair self-energy H is defined by the Dyson equation T ^ = U ^ — H. 

The dressed G and T are given in terms of the bare Go and U by Eq. and 


r = u 


l-2UGo^ 
1-UGo^ ■ 


Eliminating U yields a cubic equation for Go, which reads 

750 + 2^0 — 3^0 + 1 = 0 


in terms of the rescaled quantities 


go := Go/G 
7 ;=rG^. 

In the relevant range 0 < jqj < 2'\/3/9 the three solutions are 


( 22 ) 

(23) 

(24) 

(25) 



y97 + 4 cos 


^ 5 ( 7 ) + 27rl ^ 



I €{-1,0,1} 


(26) 


where 6{'y) = arg(—277^ — 547 — 16 — 3i7^171x74 — 277^) € (—7r,7r). The fully bold diagrammatic series 
^boid(G, T) and nboid(G, T), evaluated at the exact G(/i, G) and r(/i, G), always converge to the I = 1 
branclj^ which coincides with the exact result for |G|//i^ < (1 + \/3)/2. Above this critical interaction 
strength, the exact result is the / = — 1 branch, so that the bold series converges to a wrong result. At 
this critical interaction, the boundary of the series’ convergence disc is reached. In summary, a similar 
phenomenology occurs again. 

The consequences for the Bold Diagrammatic Monte Carlo approach [SHS] of the findings of Ref. [2] and 
of the present work is left for future study. 
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The unconventional factors 2 in Eqs. ([I5|l6|l9 ) could be removed by working with spin-dependent fi and G. 
Explicitly, we have Sboid(t^, T) = l]/t^ and nboid(t^, T) = Sboid(G', 
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